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A method to facilitate the consistent inclusion of cross-section measurements in proton parton density functions (PDFs) 
fits in NLO QCD has been developed. It allows the a posteriori variation of the renormalisation and factorisation scales 
and of the centre-of-mass energy and the inclusion of arbitrary PDFs, strong coupling in cross-section calculations. 
The theoretical uncertainties in the double differential single inclusive jet cross-section at a y/s = 7 TeV have been 
studied using this method. 



1. INTRODUCTION 

The Standard Model (SM) of particle physics, based on the theory of electroweak and strong interactions, provides a 
solid basis for measurements in high energy physics. The experiments at particle colliders LEP, HERA, TEVATRON 
and many others have proven the validity of the SM and measured its parameters. However, the SM cannot be 
considered as a complete theory, since it doesn't include gravity, the forth fundamental force, and contains too many 
free parameters, only few of them are measured to better than 1% accuracy. The new proton-proton collider LHC 
will produce a large amount of experimental data, which, due to extremely high luminosity, will have small statistical 
uncertainty on measurements for most of the production processes. However, since the center of mass energy of 
modern accelerators is far below the Planck scale, one might have to look at small deviations of experimental 
distributions from the SM predictions to find signals of new particles or interactions. The understanding of the both 
theoretical and experimental uncertainties will be the key for a discovery. 

The calculation of cross-sections at Next-to-Leading Order (NLO) in QCD involves numerical integration over the 
phase space of the final state partons in order to cancel the infra-red and collinear divergences. Since the convergence 
of such a calculation in certain parts of phase space might be slow, the integration of cross sections for observables 
at deep inelastic scattering (DIS) experiments or hadron-hadron colliders requires a time consuming generation of a 
large number of event weights in order to achieve low statistical uncertainty in theory prediction. Moreover, in the 
iterative PDF fit procedure or for the estimation of the uncertainty in the prediction due to PDFs, a s and/or scale 
^ • choice uncertainties, the full calculation needs to be repeated over and over again. 

Nevertheless, the factorisation of a cross section into long-distance PDFs and short-distance matrix elements 
(ME) allows to fill the latter into look-up tables (grids) binned in incoming partons momentum fraction and the 
characteristic scale of the production process during the time consuming evaluation and use it later on for the 



fast calculation of the theory predictions a posteriori. Several approaches exploiting this idea [l|-[3j have developed 
independently. 

The other possibility is the use of the ratio of NLO to leading order (LO) cross sections (k- factor). Such a factor, 
being calculated once, could be used to rescale the relatively fast calculation of LO cross section to mimic the NLO 
result, however, this method is only an approximation. 

This contribution describes the APPLGRID method 0] and shows its usage for the evaluation of the theory 
uncertainties in the inclusive jet cross section in proton-proton collisions at 7 TeV centre-of-mass energy. 



2. OVERVIEW OF THE METHOD 

In a typical calculation of the cross section of production process H\H2 — * Y at NLO the integration over the 
process phase space is replaced by the summation of the ME weights a over the set i = 1 ... N of kinematical 
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configurations of partons factorised in the initial state hadrons with the momentum fractions {x\, x 2 }i at the scale 



do=2^2^[ 2jT j U/bi { x u,QfJ f q2 /H 2 (x2i,Q Fi ) -tq ■ (i) 

*=1 (P.) V J 



In general there are 13 x 13 different combinations of initial state partons contributing to an observable, however, 
due to the symmetries of the ME weights, one can define a set of sub-processes and group the PDFs into "generalised" 
PDFs 

fm/m KQ 2 ) fn/H 2 {X2,Q 2 ) = (x U X 2 ,Q 2 ) , (2) 

m,n 

where index (k) denotes the sub-process. Such "generalised" PDFs are, of course, depend on the process and the 
perturbative order. For example, in the case of W— boson production one could define two sub-processes at LO and 
six at NLO 0. 

In the beginning of the evaluation of sum (QJ the three-dimensional N X1 x N X2 x Nqi look-up tables in {x\, x%, Q 2 ) 
phase space (N x and A^q2 are user defined) are created. Instead of using directly the parton momentum fraction x 
and parton resolution scale Q 2 the following variable transformation is applied 



x 



y2 



y{x)=\n- + a{l-x) ] t{Q z ) = In ( In 2_ ) . (3) 




This transformation provides equidistant binning in regions of phase space, where the PDFs are steeply falling. The 
parameter A should be chosen of the order of Aq£<tj, but need not necessarily be identical, the parameter a serves 
to increase the density of points in the large x region (A and a are user defined). 

Assumming that continious parton density function could be represented by the Lagrange interpolation 

/(*> Q 2 ) = E E / (»«. t y) (v (*) ~ w) it ] (Q 2 ) ~ ^ ( 4 ) 

i=0 j=0 

over n y and n T nodes in (y, r)— grid with sufficient precision, the sum ([I]) could be rewritten such 

da 1 

^(ttx.lfe.rj)-^ (5) 

that interpolation is performed during the event loop over the pQCD weights rather than over PDFs. Here Q 2 R — Q 2 F 
was assumed for simplicity, but this is not the fundamental limitation, since the observable scale dependence could 
be fully addressed a posteriori. 

The APPLGPJD library works as an interface, which interpolates the pQCD coefficients da^ ^ Y over {ViiV^i T ) = 
(xi, X2, Q 2 ) nodes of look-up table and fills the grid structure with one 3D-table per sub-process, perturbative order 
and observable bin. 

After the evaluation of the sum ((T|) the grid structure is written out to a ROOT file. The APPLGRID library 
provides the methods for convolution (using the Eq. [5]) of the pQCD coefficients stored in the grid with any choice 
of PDF set, a s and scales choice. Depending on the grid size and PDF set the convolution time is 1 -j- 100 ms, since 
it is already just a simple sum of real numbers and not the evaluation of divergent integrals. 

The choice of the grid architecture depends on the required accuracy, on the exact cross-section definition and on 
the available computer resources. The most critical parameter is the number of £-bins, which must be large enough 
to accommodate strong PDF variations. It has been shown [3j that the grid reproduces the original calculation with 
the precision better than 10 -3 with A^ x 20 -j- 30. 
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3. SINGLE INCLUSIVE JET CROSS SECTION 



In this section we describe the use of APPLGRID for the estimation of theory uncertainty on the example of single 
inclusive jet production of jets with radius R = 0.6 for px > 40 GeV defined by the anti-fc^ jet algorithm Q in 
several bins of jet rapidity. The cross section was calculated using NLOJET++ [5ru| with the renormalisation and 
factorisation scales equal to the pr of the leading jet. 

Figure [1] shows the relative theory uncertainty in the cross section. Both the total uncertainty (orange band) 
and all its individual components are shown. The scale uncertainties (right side hatches band) are defined as the 
envelope of the independent two scales varying by a factor two up and down of the default choice. The effect of the 
uncertainty on the a s determination (yellow band) on the observable is estimated by calculating the cross section 
using strong coupling values within the uncertainty range and using PDF sets fitted using these values 12j . The PDF 
uncertainty (left side hatches band) is calculated using the uncertainties PDF sets as recommended by the authors 
If! 11 , 2, 14 [ . Finally, the individual components are added in quadrature in the total uncertainty. One can see that 
both CTEQ6.6 and NNPDF2.0 give similar uncertainty for the cross section in the central and the forward region 
for jets with px < 100 GeV. However, the use of NNPDF2.0 allows to shrink the total theory uncertainty to 
compared to CTEQ6.6 estimation of ±11% for high pr central rapidity jets. In the forward region the difference 
between the CTEQ6.6 and the NNPDF2.0 prediction for high pt jets is even more pronounced : up to for 
CTEQ6.6 compared to ±? 2 % for NNPDF2.0. 

Figure [2] shows the comparison of different PDF sets predictions in several bins of rapidity (up to \y\ < 2.8). Also 
the predictions of all the PDF sets agree between each other within the error bands, the shape of pr dependence is 
quite different for different PDF sets. All the PDFs agree within 4% at low (pr < 100 GeV) transverse momentum, 
but NNPDF2.0 and MSTW2008 predict 5-7% more jets at intermediate p r 's, while HERAPDF1.0 estimates 5% 
less cross section in this p T region. At high jet transverse momentum HERAPDF1.0, NNPDF2.0 and MSTW2008 
predict, depending of rapidity region, a ~ 15 — 30% smaller jet rate compared to CTEQ6.6. The NNPDF2.0 and 
MSTW2008 provide very similar pr dependence over all the rapidity bins, they agree to CTEQ6.6 at low pr, predict 
larger cross section than CTEQ6.6 at intermediate pr values and lower one for high pt jets. The cross section 
calculated with HERAPDF1.0 is always smaller compared to one with CTEQ6.6, the diference is small at low pt's 
and it increases towards higher transverse momentum. 



4. SUMMARY 



APPLGRID is an open source software project written in C++ which provides a framework to create production 
process dependent look-up tables of perturbative coefficients for a posteriori calculations of observable with any given 
PDF set, choice of the scales and strong coupling. The method is based on the Lagrange interpolation of PDFs and 
the use of symmetry in hard scattering allowing to group weights into sub-process contributions. It offers a robust, 
transparent and uniform way to study the theoretical uncertainties for various QCD and electroweak processes via the 
interfaces to jet production cross section calculator NLOJET++ 0-0] and MCFM @,@|. It provides the possibility to 
look at contributions of different sub-processes to the observable. APPLGRID is a very efficient substitution for the 
use of k-factors in PDF fits to the jet and electroweak cross sections data. The software offer additional functionality, 
such as a posteriori centre-of-mass energy rescaling, arbitrary variation of renormalisation and factorisation scales 
and the interface to use fastNLO [l| grids for DIS and hadron-hadron collisions. 

The use of APPLGRID for a posteriori evaluation of uncertainties due to renormalisation and factorisation scale 
variations, strong coupling measurements and PDFs uncertainty sets has been demonstrated on the example of single 
inclusive jet cross section. The APPLGRID has already been used as the standard tool by the ATLAS Collaboration 
for the interpretation of the jet measurements 15| • 
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Figure 1: Theory uncertainties ifor the double differential single inclusive cross section as a function of jet transverse momen- 
tum. The two upper (bottom) figures compare uncertainties calculated with CTEQ6.6 and NNPDF2.0 PDF sets in central 
(forward) jet production. 
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Figure 2: Comparison of MSTW2008 Q (blue circles), HEPAPDF1.0 [IH (black triagnles), NNPDF2.0 (violet squares) 
and CTEQ6.6 12] (yellow band) predictions of single inclusive jet cross section and the total theory uncertainty uncertainty 
normalised to the CTEQ6.6 central prediction as a function of jet transverse momentum . 
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